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ABSTRACT 


This report discusses the design and implementation of 
a fixed-based spin simulator and the results derived from 
conducting preliminary spin tests on the simulator. 

The central piece of equipment in the simulator was a 
hybrid computer in which the analog computer solved the 
equations of motion while the digital computer performed the 
tasks of program control and aerodynamic data storage. The 
visual display consisted of a computer-drawn picture on a 
graphics terminal, while pilot control was obtained by use 
Sema simulated cockpit situated in front of the graphics 
terminal. : 

Results showed that the simulator displayed excellent 
dynamic response characteristics and provided sufficient 
visual cues to perform meaningful spin tests. 

This project was a continuation of previous work and has 
shown that the design and construction of this simulator has 
been an excellent research tool and source for further study 


in the field of control systems and aircraft dynamics. 


aia 
. yo A i- 7 cs : . 
2 _ _ 


> os i 
~ - . 
- - ~~ * = - 4? 
ue ee ae 
a a en 


.- 





TABLE OF TCON TENTS 


Tf INTRODUCTION ---~--~------------------------------- 9 
II. THE SIMULATOR ----------------------------------- 12 
III. SIX DEGREES OF FREEDOM AIRFRAME EQUATIONS ------- 15 
IV. GRAPHICS ---------------------------------------- 20 

A. GRAPHICS PRESENTATION ----------------------- 20 

B. GRAPHICS PROCESSING ------------------------- 21 
Vv. RESULTS OF SPIN TESTS --------------------------- 29 
VI. CONCLUSIONS ------------------------------------- 36 
APPENDIX A - THE DIGITAL PROGRAM ---------------------- 37 
APPENDIX B - THE ANALOG PROGRAM ----------------------- 69 
APPENDIX C - PROGRAM VARIABLES ----------------------- - 79 
APPENDIX D - DATA DECK PREPARATION -------------------- 88 
APPENDIX E - OPERATING MANUAL ------------------------- 92 
APPENDIX F - SAMPLE OUTPUT ---------------------------- 96 
LIST OF REFERENCES ------------------------------------ 135 
INITIAL DISTRIBUTION LIST ----------------------------- 136 
FORM DD 1473 ------------------------------------------ Sy 








“ 


piSi OF Ee iGures 


Simulator Control Loop ---------------3- 9-9 --------- 10 
Aircraft Coordinate Axes -------------------------- 16 
Earth Grid Reference Lines ------------------------ Ze 
Initial Graphics Display -------<<<---.<<<---------- Zo 
Graphic Coordinate System ------------------------- 24 
H-Matrix ------------------------------------------ 29 
Euler Angle Rotation ----------------------------- - 26 
Results of Spin Test ------------------------------ il 
Representative Spin Results ----------------------- 55 





biotsOn sovMiBbOLS 


WiIMG “Spal ase es 
Rolling moment coefficient 


Lote nOmeecOnetolecthectlvyeness derivative 


ROiMng WOMent COSLiEbcLent due to rudder deflection 


Waliping In Toll derivative 

Rolling moment coefficient due to yawing 
Pitching moment coefficient 

Elevator control effectiveness derivative 
Pitch damping derivative 

Yawing moment coefficient 

Aileron yaw derivative 

Rudder control effectiveness derivative 
Yawing cross derivative 

Damping in yaw derivative 
Bongteucinal=-seoree cockhlerent 


Longiteucmnal—fhoree cocikiicient dcue-to elevator 
deflection 


Longitudinal-force coefficient due to pitching 
Side-force coetflrcient 

Side-force coefficient due to aileron deflection 
Side-force coefficient due to rudder deflection 
Side-force coefficient due to rolling 

Side-force coefficient due to yawing 


Vereical—-tforce cocrtficient 


OO} 


Vererecal—-force coeriicient due to elevator deflection 


Vertical-force coefficient due to pitching 
Wien Chorcus ft, 

Aerodynamic force in X body axis direction, lb. 
RO HOReCe wie wis tabtitty axis direction, 1b. 
Hotelenorcecwmn « wimd axils direction, lb. 
herOdyMamee £OrCce im x body axis direction, Ib. 
iota worce in Yo stability axis direction, ib. 
Hoe weenoOrece ime wind /axis direction, lb. 
Aewedynamic ELorce in 4 body axis direction, 1b. 
Total force in 5 stability axis direction, lb. 
Ute rOrce Ane 4 Wwithe ia<1S idaerection, ib. 
Acceleration due to gravity, 32.2 ft/sec” 


Moment of inertia about X body axis, slug-£t* 


Gross, PEOaUUuct: Of inertia, slug-ft* 


Moment of inertia about Y body axis, atc eae 


Moment of inertia about Z body axis, slug-ft* 
Rolling moment about body axis, ft-lb 
Puteching moment about body axis, f£t-ib 

Mass of aircraft, slug 

Yawing moment about axis, ft-lb 

Rolling rate about body axis, rad/sec 
Rolling rate about stability axis, rad/sec 
Normalized rolling rate about body axis 
Pitching rate about body axis, rad/sec 


Pitching rate about stability axis, rad/sec 


Normalized pitching rate about body axis 





q Free stream dynamic pressure, Ney 


R Yawing rate about body axis, rad/sec 
R. Yawing rate about stability axis, rad/sec 
R Normalized yawing rate about body axis 
S Aerodynamic reference area, feo 
S : 
x 
S Y 
y 
S Z 
Z, 
ah Tieust oom aaLrCrart, “ib. 
V Velocity of aircraft, ft/sec 
X “anerpelal, CoOorLaginate Of Sirecrart, ft. 
a MbnnCtiire ime COOLdLNakeeOr a1rerart.. Lt. 
Z, BEeInercr lal iCOOEGI nace, OF clrerart..» Lt. 
a Angle of attack, deg. 
B Angle of sideslip, deg. 
6a Aileron deflection angle, deg. 
Onde Elevator deflection angle, deg. 
One Rudder deflection angle, deg. 
) Euler roli angle, deg. 
Uy Fuler yaw angle, deg. 
6 Pile tapiteh «ang le, Geq. 
i) SHON she 








ACKNOWLEDGEMENTS 


The author wishes to express his sincere appreciation 
to: Lieutenant (junior grade) Michael Redlin, Assistant 
Professor of Aeronautics, for his invaluable guidance, time 
and assistance during the course of this project; Robert 
Limes, William Thomas and Al Wong, of the Department of 
Eileeccical Engineering Computer Laboratory, for their 


assistance in the .running of the program. 





PLN RODUCTION 


Men puioOceeOl this project was a COntinuation study of 
the work done by J. H. Kahrs [Ref. 1]. Kahrs made the 
initial design and construction of a fixed-based variable- 
stability simulator with the task of landing on a carrier or 
runway. The intent here was to make the necessary modifi- 
cations and adjustments to this simulator in order to 
develop a simulator capable of evaluating the spin charac- 
betristics Of various aircraft. 

The same basic equipment, l1.e. the hybrid computer, 
graphics terminal, and the cockpit simulator, was used as 
constructed by Kahrs. The hybrid's analog computer was used 
BOmethne SOlution Of the equations of motion of the aircratit, 
while the hybrid's digital computer was used for overall 
program control, graphics generation, data reduction and 
storage. The visual display was created by a graphics 
processor and terminal. The cockpit simulator, situated in 
meeont Of the graphics display, consisted of a chair with 
Various control levers and buttons whose outputs were tied 
directly to the analog computer, thus giving the operator a 
@grect means of controlling the simulator. This combination 
of computers and Somerol linkages comprised the basic hard- 
ware of the simulator and formed a complete control loop as 


diagrammed in Figure l. 
















xDS-9 300 
DIGITAL 
COMPUTER 


Er—s000 


ANALOG 
COME ara 





ADAGE AGT-10 
DER =Z 
GRAPHICS 

PROCESSOR 





COCKPIT 


Prguices ta Simulator Control hoop. 


Me 





The simulator, in the above format, provided the neces- 
sary and required operating and handling characteristics. 
These characteristics, as originally designed by Kahrs, 


maclmuaded (1) simplicity of setup and operation, (2) ability 


memenenge alrcrart being simulated, (3) capability of simu- 
imaeing non-linear aerodynamic data, (4) automatic scaling of 
the analog computer solution, (5) fixed-base inside-out 


display, and (6) ability to expand or change simulator task. 

| The two most important characteristics desired in this 
project were the ability to expand and change simulator task 
and the ability to change aircraft being simulated. Major 
revisions in the internal programming were needed to account 
for the large differences in the aircraft's flight character- 
istics between the landing phase and the spin phase. These 
differences included the addition of the non-linearized 
equations of motion that are needed in the solution of the 
spin problem. The ultimate goal was to verify that this 
simulator was capable of providing a realistic simulation of 
a spin from a pilot's point of view and at the same time 


obtain meaningful data for further study. 
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iin sli oO LMULATOR 


The simulator consisted of three basic pieces of equip- 
ment; the hybrid computer, the graphics terminal, and the 
cockpit simulator. The hybrid computer was used due to the 
outstanding qualities obtained by use of the digital computer 
in conjunction with the analog computer. Due to the rapid 
computational time of the digital computer the bulk of the 
work load was assigned to it. This included the resolution 
of the forces and moments, build-up of graphics data, and 
the non-linearities of the aerodynamic coefficients. The 
analog computer, on the other hand, was used for integration 
@metne equations of motion and interfacing. This also 
provided the most rapid means of obtaining a desired output. 
The important feature involved throughout was obtaining the 
desired output data by the fastest means possible. fThis 
assured that the simulator had an excellent dynamic response 
and led to a display with as little flicker as possible and 
Overall grapnic quality. 

The graphics terminal consisted of a graphics processor 
and visual display unit (CRT). Construction of the display 
data was accomplished by the digital computer. These data 
were sent to the graphics processor where a computer-drawn 
picture was formed on the CRT. As new data were compiled by 
the digital computer, the graphics processor would update 
and refresh the visual display, thus providing the pilot 


Wiech, a CONLINUOUS DIcture. 
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The cockpit simulator was a chair, with attached 
Soenmterols, placed in #£ront of the CRT. A Gemini control 
stick was attached to the right arm and provided all three 
Genero! movements: yaw, pitch, and roll. Two programs 
control buttons were provided. One control button was 
located on the top of the stick, while the second control 
button was located on a throttle plate attached to the left 
emnalong with a small throttle control. All outputs of the 
Geckpit were connected directly to the analog computer for 
processing, thus providing the pilot with the necessary 
controls to run the problem from his station. A complete 
operating manual 1s contained in Appendix E. 

The simulator was designed to accept the following input 
@aeaein the form of punched cards: aerodynamic coefficients; 
Mmepal conditions (initial position, altitude, velocity, 
etc.); aircraft constants (weight, moments of inertia, wing 
span, etc.); scale factors based on maximum expected range 
of the problem's variables; and an earth reference grid. 
Appendix F outlines the preparation of this data deck. 

The simulator then automatically processes these data, 
sets the potentiometers of the analog computer, and sets the 
initial conditions on the integrators. Control is then 
transferred to the cockpit where the pilot has three options 
available: 

Po hOommlvworie= cimulator with the further, option, oF 

aborting the run at any time and returning to fly 


again with the same aircraft parameters. 


Ie. 





2. To input a new data deck to change the simulated 
Grueckart Sspamatoters Or to change aircraft» being 
Simulated. 

3. To stop the complete problem. 

At the completion of each run an analysis of the run is 
briefly displayed on the CRT after which instructions for 
further program control are displayed. The analysis of the 
run consists of the following items: 

ime Pinel angle of attack 

2. Final pitch angle 

Sin eileen seeel it angle 

4, Final yaw rate 

Bee final velocity 


(reine! altitude 
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Tit. SIX DEGREE OF FREEDOM ATRFRAME EQUATIONS 


Three different coordinate systems--body, wind, and 
stability axes--were available for problem solution. It was 
Gemcluded, on the basis of work done by R. M. Howe [Ref. 2] 
mie. H. Kahrs [Ref. 1], that the best choice of the coor- 
dinate systems was a combination of the wind axes system and 
the body axes system. The translational equations of motion 
would be based on the wind axes system while the rotational 
equations were based on the body axes system. The three 
coordinate systems are shown in Figure 2. The stability and 
body axes system differ by the aircraft's angle of attack, 
while the wind and stability axes system differ by the side- 
Slvp angle. 

With a knowledge of the three angular rates of motion of 
eiemabrecratt (P,0,R) about the body axes, normalized values 


feorer, ©, and R are obtained. 


- Pb 
p= oa (1) 
mas (lo) 
- Rb 
ae hy 


The angular rates about the stability axes are determined as 


monet lons 1or PP, O,--and Rk. (Eqn. 4-6) . 
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P cosa + R sina (4) 
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Figure 2. Aircraft Coordinate Axes. 
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Re = P sing + R cosa (6) 


The aerodynamic forces expressed in the body axes can 


now be calculated. 
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The body axes forces (Eqn. 7-9) can then be resolved into 
forces along the stability axes. Thrust (assumed to be in 


the X-direction only) and the force of gravity are also 





inserted. 
F F 
De ae 
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2aero 
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The final transformation of forces from the stability 


axes to wind axes can then be made (Eqn. 13-15). 
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The translational equations of motion can now be found 


by solving for the derivatives. 
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ry ar 
V=—-— (16) 
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The Euler angular rates are also solved now (Eqn. 19-21). 
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Next the moment equations about the body axes are solved 


fFOorNas  Snown below. 
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The rotational equations of motion can now be written 


Migm. 25-27). 





oC ee I (R+PQ) 
yy S YY 27, xz L 
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In the development of the velocities in the inertial 
frame, a small angle approximation was not used due to the 
fact that large angles are frequently encountered during the 


Spoeaeproblem. The resulting equations are as follows: 
S - v|cosucostcostcos) + sing (-cos¢siny + sin¢sintcosv) 


+ sinacos8 (singsiny + cos¢sinécosy) (28) 


§ = v|cosacosscostsinv + sing (costcost+singsintsiny) 


4: sinasing (singcosy + cos¢sinésinv)| (29) 


S, = v|-cosacosBsind + sing (sinécos®) 


= Sinacoss (costeos#) | (30) 


Integration of the various derivatives derived in the 


above equations yields the desired state variables. 
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iV. eGhAPHLCS 


A. GRAPHICS PRESENTATION 

The need for a good graphics presentation was evident 
due to the complete lack of any physical and audible cues to 
the pilot. Therefore a graphics display had to be generated 
which would give a representation of the earth's surface and 
provide the pilot with enough visual cues to enable him to 
orient himself throughout the problem. The essential visual 
cues needed for the pilot's orientation were indications of 
roll, pitch and yaw. A horizon was included for the orien- 
jeeeron in pitch and roll. Due to the continuous yawing 
motion of an aircraft during a spin, it was vital to have a 
graphics picture that would accurately display this motion. 
After investigation of several possibilities, it was con- 
cluded that a grid system, attached to the earth's surface, 
composed of a one mile per side squares would provide the 
necessary requirements. The entire grid system was limited 
to six squares per side to avoid excessive computational 
time. By extending the grid lines at all four cardinal 
headings, the visual representation of the yawinog motion of 
the aircraft was further enhanced. This improvement was 
accomplished by the fact that during a continuous yaw the 
pilot only sees a series of lines every 90 degrees. He is 
thus able to positively determine the yaw direction as well 


as a relative measure of the yaw rate. These extensions of 
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the grid lines also assist in the orientation of pitch and 
molb, Ihis grid system is shown in Figure 3. 

To provide the pilot with another source of visual cues, 
a Single instrument was constructed. This instrument pro- 
waraed the pilot with two sources of information. The first 
was a "needle" providing the pilot with the aircraft's angle 
@mebank and the second was a “ball” providing the pilot with 
a Measure of the aircraft's sideslip angle. This instrument 
was not the normal "needle-ball" found on all Naval aircraft 
but was a combination of an attitude gyro and turn and bank 
Mmmeeecator,. It was felt that the angle of bank information 
taken from the attitude gyro combined with the sideslip 
angle information taken from the turn and bank indicator 
provided the pilot with the best information in the limited 
Space available. 

In the final display, the picture that the pilot Saw was 
Pestm1eted in direction to the X-axis of the aircraft and to 
a square with a field of view angular limit of plus or minus 
18.5 degrees. Within this square or window lies the horizon 
and portions of the grid reference system. As the aircraft 
maneuvers in space the horizon and grid reference system 
move dynamically within this window creating the sensation 
of flight. The complete initial display is pictured in 


Figure 4, 


Doone bear ROCESS ING 
The processing of the grid reference system from a three- 


dimensional object to a two-dimensional plane capable of 
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Figure 3. Earth Grid Reference Lines. 
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Figure 4. Initial Graphics Display. 
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being displayed on a screen proved to be the most complicated 
portion of the display. The same system developed by 

L. G. Roberts [Ref. 3] and R. B. Desens [Ref. 4] and used by 
J. H. Kahrs [Ref. 1] was used here. It basically consists 
@amene CONSErUCEION Of a single transformation matrix for 

all display points. Two coordinate systems were used in the 
development of this matrix, that of the earth reference 
system (object) and that of the aircraft (viewing plane), 

and are shown in Figure 5. 


fhe teanstormatioOn Matrix (H-Matrix) is the product of 


five matrices: rotation, translation, perspective, scale, 

and another translation. The H-Matrix is further reduced 
Viewing X 
Plane P 





Figure 5. Graphic Coordinate System. 


24 





by Desens to the product of three matrices as shown in 





Eeecuire 6. 
AG By JS i o/c 0S 7E 
By Ba B, 0 it 0 0 
he = 
Cy C. C4 0 0 it 0 
Am Yn Zin 0 Yo Zo S. 
Object Axis Viewing Plane Perspective 
Rotation Orientation 
Scale 
Viewing Plane 
Translation Offset 
ba gGume 6ee) H=Matri x. 


Whe f£127St matrix consists of the direction cosines of 
the Euler angle rotation of the object axes and the trans- 
H=etom distances between the two axes system. The direction 
cosines and the order of rotation are shown in Figure 7. 

Since the object axes in this problem was a fixed earth 
reference system, the direction cosines are all zero except 
for the diagonal terms which have a value of one. The 
matrix then is simply reduced to a viewing plane translation. 


Mae translation distances; X Y. and Zin are derived as 


faire fe 
Pol LOws : 
Xm = Aopject ~ *view plane 
Yep 7 ‘object 7 view plane 
7, —_ 


‘oe “object 7 


Tan 
view plane 
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A. = cosfcosa 


il 
A, = cosfsina 
A. = - sing 
By = - cosysina + sinysingcosa 
Ba = cosacosy + sinysinfsina 
B. = sinycos®f 
Z 
Cy = sinySina + cosysinfcosa 
O 

Cy Seo tl) COSA eos y si Nhsing Y 

=P Bp 
C3 cosBCcosy ‘ 


OrdertoOrmene cation. Cre ,y 
Figure 7. Euler Angle Rotation. 
The second matrix 1s composed of the direction cosines 


of the Euler angles of the aircraft or viewing plane axes 


and are listed below. 


AA, = cosycosé 

AA, = sintcosd - coswsin@sind 
AA. = sintsind + costsin@cosd 
BB, = - cos@sinwt 

BB, = costcosd - sin@sinvsind 
BB, = costsind - sintsin€cosd¢ 
CC, = - sing 
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CC. = - cos8sind 
CC. = cos8coso 


Ordervor Gotat lon 


Wivrawo ere (ottcn) , droll) 


The third matrix is composed of the offset option, Zo 
and MOE the scale factor (S) and the focal length(F). The 
offset option was not used but its usage is explained in 
Ref. 4. The scale factor was set at 1/2 and the focal 
femgen, the distance between viewer and the viewing plane, 
was set at 1.5 ft. 

To process the grid reference system or object, a point 
(e752) 15 taken and converted to homogeneous coordinates by 
the addition of a scale factor (W). This yields new coor- 
peices (X',Y',Z2*‘,W) where X' = WX, Y' = WY, Z' = WZ. Post 
multiplying these coordinates by the H-Matrix yields 
Can 72 ,W'). The display coordinates can now be found by 


Gaviding by W'. 


ue 


It 


y" /W' 
I = z" /W! 


The X-coordinate was not used in this display since it 
is only an indication of Ccepth. The final displav coordin- 
ates are obtained by passing each point of the grid reference 
Syoecem through the M-Matrix. 

The horizon was constructed by drawing a line parallel 
to the Y-object axis and offset a distance in the negative 


M-Garectron. The Z-cooradinate. was set to. zero. itn -doing sce 
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the assumption was made that the earth was flat and there- 
meee did not account for the fact that when flying at an 
altitude a depression angle 1s created due to the curvature 
of the earth. This assumption was valid since the horizon 
was only used for a relative measure of orientation for the 
pilot and not as a measure of the altitude of the aircraft. 
The coordinates of the horizon are passed through a modified 
Viewing plane orientation matrix to account for pitch and 
meme Of the aircraft. The yaw angle of the aircraft is set 
to zero since it does not have any apparent visual influence 
on the movement of the horizon. This results in the follow- 


migematrix. 


AA, = cosd 


AA, = - sin8sind 


AA. = sin@cosd 


BB, = 0 

BB. = cos$¢ 

BB, = sind 

CC, ns? eels 
CC. - cos@sind 


CC 3 = cos8cosd 
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Vee obo Or vor EN TROLS 


Several spin tests were conducted in order to compare 
these results with the results of other simulators using the 
Same aircraft parameters and aerodynamic coefficients. The 
aircraft parameters and aerodynamic coefficients used in all 
tests were eer eee W. P. Gilbert [Ref. 5] and are 
representative of a variable-sweep fighter aircraft. A wing- 
Sweep angle of 16° was the only configuration considered 
during the tests. 

In all cases the aircraft was initially positioned at an 
altitude of 30,000 feet and an airspeed of 622 ft/sec. The 
spin was entered by reducing the throttle to the idle posi- 
tion followed by movement of the horizontal stabilator to 
miyemwrull trailing edge up, rudder to full trailing edge left, 
and the ailerons to full right wing down. The controls were 
Maintained in these positions until the spin had fully 
developed when they were then released to the neutral posi- 
tion. The resulting spin was limited in altitude to 10,000 
feet or until the spin development was such that the limits 
of the computer were exceeded. At no time was an attempt 
made to effect a recovery. 

The resulting aircraft motion was a stall exhibiting 
Mie tsablestligit characteristics in lateral control. ~A._right 
St letteroll of 360 at high angles of attack then paeoes 


a@eer wiich the abrerart entered ‘a fast flat-spin. During 
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the spin, the average angle of attack was 80° with a yaw 
rate of about -133 deg/sec. The pitch angle and roll angle 
remained relatively constant throughout the spin at approx- 
Meee ly —15° and 0° respectively. At the end of 60 Seconds, 
the aircraft had completed 9.5 turns and lost about 7000 
feet of altitude. The final aircraft velocity was 262 
ft/sec. A complete time history of the results of one of 
the spin tests is shown in Figure 8. A brief comparison 
with the results obtained by Gilbert is shown below. A 


complete time history of Gilbert's results is shown in 


Figure 9. 
Simulator Gilbert 
1. Average a 80° Cor 
2. Yaw Rate -133 deg/sec -160 deg/sec 
3. Average @ Sl Sve ile 
4. Average $¢ Or On 
Se Number Turns el: Oo 
6. Altitude Lost 7000 feet 8000 feet 
ten Length Of Spin 60 sec 40 sec 





ieee 


b 5) 


Zk 











hi, 
1 
nti 
{ 
ia : 
! 
t 





it! 
“1 leit i! 
Oe he | 
a pe 
= ine ons 
/ Se any 1. (| 
reg ER RPERA ua ce , Wn, ° 


oO oOo OF BW 








50 60 


40 


30 
Time, 


20 


10 


sec 























20 40 50 60 
Time, 


20 


10 


SecC 


\ } 
HN 
- 


ane 
VAN 


\ 


‘ + 
4 
} 
i} 
iN = 
==. a 
+ 
t Me } 
f=-— $ ~- 





=EEEECENAN 





AN 
30 | 
Time, 


mia 
a 
ETN 


a 





tt 
Fa 
- 
| | 
| 
et i ik malic | \! 
is { ty | 
th ill 





60 


50 


40 


20 


10 


Sec 


Results of Spin Test. 


Figure 8. 


gl 





60 



















































































7 a 
ney 
|e 
tin i i!| Lt) 
eT EY 
5 | : | » Teeny Pet 7 sw 
i EE A : 
o EEE He! 3 
49) a! WY 
) tifa 
Toh pos ) ~ 
Sa ita! ! 8 e 
4) aiid | 
= To | 
a HH 
A ae a 
/| 
SD a 
eA N 
PR 
Te AT eh 
REE EAN a 
LH A A 4 
DEF A NR 
© 











30 
Lime. Sec 
(Cone ilnued).. 


20 
32 


10 
Brgure 8. 





fra ie ER 
ee ee EE 
SSE 
SEO ein aN 
re 








30 40 50 60 


Time, 


20 


10 


sec 














‘ ' 
\¢ ' 
1d 





30 40 50 60 
Time, 


20 


sec 








PCE eer 
Sone Oat 





188 


SseC 


Time, 


(Continued) . 


Figure 8. 


oie) 











30 40 50 60 
Time, 


20 


10 


sec 



















































































Slee 


Time, 


(Conernueda)- 


Figure 8. 


34 





r-) 
a 
3 


8 


TTT 
= 
o 8 
: 


B, deg 


2 86 8 £8 BS 8 8 


0 —— 100 — 


150 deg/sec . 











q 
OQ ———— — ——— - ---- ee - 
deg/sec . 
ay. 4 - 
g units 
0 t— \ Ms 
Right §$ 
19) a \ 
' : 
Turns of Vp. 2 400 
mf sec 100 -— 
= ce — 
left 5S 1 —-|200 
59 — 
be o— 0 





Time, sec Time, se 


Figure 9. Representative Spin Results. 


35 





Vi eeONncEhuS TONS 


As a result of conducting a series of spin tests on this 
Simulator the following conclusions have been reached: 

1. The simulator provided excellent visual cues to 
enable experienced pilots to orient themselves 
throughout the entire problem. 

2. The simulator exhibited good dynamic response with a 
computational loop time of 80 milliseconds and solu- 
tion update rate of approximately 16 samples per 
second. 

3. The simulator was capable of providing meaningful 
Spin data from which comparisons and further study 
could be made. 

4. One minor drawback to the simulator was the lack of 


any physical cues, i.e. "g" forces. These forces can 
be significant in an actual spin and severely res-~ 
trict the pilot's movement. Lack of these forces 
has led to pilots over~-controlling the simulator 
and thereby risking the possibility of introducing 
erroneous information or exceeding the limits of the 
computer. 

5S. ihe versatility of the simulator, in the ease with 
which the aircraft or aircraft parameters may be 
changed, has made it a valuable research tool and 


Source tor further study an/the field of flight 


dynamics. 
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APPENDIX A 


THE DIGITAL PROGRAM 


This appendix contains a listing of the computer program 
written in the FORTRAN language. It consists of a main 
program and three subprograms: the data reduction, the 
computation loop, and the A/D-D/A. The A/D-D/A subroutine 
Mieritten in METASYMBOL, the assembly language for the 


XDS-9300, and was prepared by the computer laboratory staff. 
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PROGRAM VARIABLES 


The following is an alphabetical listing of the program 


variables used in the computer program. 


A Array of angles, in order, alpha, beta, psi, theta 
and phi 

Al10 ideal Conditions on A(l), the vanglerot sattack 

AA Array of digital to analog variables 

AD Array of analog to digital variables 


ADOTN Array of the negatives of the time derivatives of 
the angles inA 


ADOTO Old derivatives of angles in A used for prediction 
ALI Aerodynamic rolling moment about the body axis 


ALPHA Real argument based on angle of attack used for 
coefficient lookup 


AMI Aerodynamic pitching moment about the body axis 

ANI Aerodynamic yawing moment about the body axis 

ARG Argument used during sine and cosine lookup (type 
real) 

B Wing span 

BZ One half the wing span B 

BETA Real argument based on sideslip angle used in 


coefficient lookup 


BX A measure of sideslip angle used in the generation 
Gite "bali. 


G Array of aerodynamic coefficients for a given angle 
of attack and sideslip angle 


CA Reway Of cosines Of the angle array iA 


719 





CB 
CBZ 
Gen (n) 
cHoeyC 


COEFA 


COEF AB 


Gor (n ) 


CON1 


FPYAM 


FYoM 


FYWM 


Mean aerodynamic chord 

One half the mean aerodynamic chord CB 

Scalar equivalents of the array CA 

Choice of program options offered by main program 


Array of coefficients dependent on angle of attack 
enly, | 


Array of coefficients dependent on angle of attack 
and sideslip angle 


Summation variables used in the interpolation 
process of coefficient lookup 


Program constant equal to the number of degrees in 
a radian 


Avleron deflection 


Degrees of angle of attack for which COEFA and 
COEFAB are tabulated 


Degrees of sideslip for which COEFAB are tabulated 
Flying tail deflection 
Initial flying tail deflection setting 


Scale factor divisor used to divide out the scale 
factor of the homogeneous transformation coordinates 


Rudder deflection 
Focal length of viewing image 


Aerodynamic force in the X-direction divided by 
the mass of the aircraft 


Force in X-direction of the stability axis divided 
by the mass of the aircraft 


BRerodynamic Force in the Y-direcerionm my ided voy 
the mass of the aircraft 


Force in the Y-direction of the stability axis 
divided by the mass of the aircraft 


Force in the Y-direction of the wind axis divided 
by the mass of the aircraft 
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FZAM 


FZWM 


HeNSCT 


HM 
a 
LA 


IAGN 


IALPHA 
IARG 
Ips 
iene A 
meno YC 


BE 


DEV 


1iJOENS 


IER 
1h eps 


oe 7 


Ho1@) 


Eo rR. 


PE xnn 


Aerodynamic force in the Z-direction divided by 
the mass of the aircraft 


Force IMetne G-adirecticon of Ene wind alc main ded 
by the mass of the aircraft 


Acceleration due to gravity 
Scaled acceleration due to gravity 


HOmizontal intersection of a line, uSedeim tie 


software window 

Viewing plane orientation matrix forthe horizon 
Integer variable usually a counter 

Do-loop counter based on angle of attack 


A counter to check the number of times a point has 
gone through the software window 


Integer conversion of ALPHA 


Integer conversion of ARG 


Do-loop counter based on the sideslip angle 


Integer conversion of BETA 


Integer conversion of CHOYC 


Graphic array containing the dynamic portion of 
the display 
Graphic device number (1 or 2) 


Grape block directory for the drapiteomeug teal 
processor 


Error flag returned by graphic subroutines 
Integer counter equal to I plus one 


Integer scale factor based on thrust divided by 
mass 


Array/graphics block of fixed data, basically the 
square or window 


Text directory for the graphics digital processor 


Line of text that is sent from the digital computer 
SOREN rapite alo tale processor 


81 





4 
7 
; 

: 








QUE 


Operating values of start/end points (X-coordinate) 
used in software window 


Equivalent of XSTART, used in software window 
Equivalent of XEND, used in software window 


Operating values of start/end points (Y-coordinate) 
used in software window 


Equivalent of YSTART, used in software window 
Equivalent of YEND, used in software window 
Be-Veop counter 

DBo-loOoOp counter 

Atbay sof integers used in outputing coctticrents 


Array of integers used in outputing analog pot 
settings 


Clock variable used in loop timing 


Number of analog to digital conversion variables/ 
truck lines 


Number of digital to analog conversion variables/ 
truck lines 


©ctal value used to null™out a lone ot text 
Variable used for program flow and eanon 
Boewlar velocity Of the alrcralt about -ai- 
P normalized 

The negative of the time derivative of P 
Peay “Or DOL Settings: Lox @ere sanalog 

P about the stability axis 

Anoular velocity Of the alreratt about tne Y arms 
Q normalized 

The negative of the time derivative of Q 
Orasout the Stability axis 


Dynamic pressure 
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Sie & 


SF4 


SiS 


SF6 


SE/ 


SF8 


Angular velocity about the Z-axis 

R normalized 

The negative of the time derivative of R 

Density of air 

Moment of inertia of the aircraft about the X-axis 
Product of inertia about the X and Z-axis 

Moment of inertia of the aircraft about the Y-axis 
Moment of inertia of the aircraft about the Z-axis 
Mass Of the alrerart 

R about the stability axis 

Graphic transformation matrix or array 

Wing area of the aircraft 

Array of sines of the angle array A 


Array of time derivatives (velocities) of the 
distances in all three inertial directions 


Product of the scale factor and focal length 


Scale factor based on maximum velocity in 
X-direction 


Scale factor based on maximum velocity in 
Mediusection 


Scale factor based on maximum velocity in 
Z-direction 


Scale factor based on maximum angular velocity in 
X-direction 


Scale factor based on maximum angular velocity in 
Y-direction 


Scale factor based on the maximum angular velocity 
in the Z-direction 


Scale factor based on thrust divided by mass of the 
alGenra Lt 


Scale factor based on maximum angle of attack 


yes 





oro 

SF10 
em 1. 
Srl 2 


Sral3 


Sr 14 


=) did lige. 


1.6 


SE 1l?7 


SEs 


Sr 19 
E20 
ee 1 
Sh 22 
Sb 3 
SF24 
5) vate. 
SE26 
S27 
S28 
Sra, 
SE 30 
Sod 


oe 32 


Scale 


Scale 


Scale 


Scale 


Scale 


factor 


factor 


bactreim: 


factor 


PAC tOn 


X-direction 


Scale 


EQeCCor 


Y-direction 


Scale 


Faccor 


Z-direction 


based 
based 
based 
based 


based 


based 


based 


Scale factor based 
deflection angle 


Scale 


FACEOL 


tion angle 


Scale 


factor 


tion angle 


Scale 
Scale 
Ratio 
Ratio 
Ratio 
Ratio 
Ratio 


Ratio 


if aAGt Ole 


factor 


based 


based 


based 


based 


Ge 
et 
Gil 
Gr 
of 


of 


Sl 6 


Sri 


SE il 


Sra 


ord 


ey es, 


SF4 squared 


SF5 squared 


SF6 squared 


to 


to 


to 


to 


EO 


1, 


SF4 


SF4 


Ses, 


SEG 


Sr 


oi 6 


OM 


on 


On 


on 


On 


(e168! 


on 


on 


on 


on 


on 


on 


the 


the 


the 


the 


the 


the 


the 


the 


the 


the 


the 


sig 


maximum 
maximum 
maximum 
maximum 


maximum 


maximum 


maximum 


maximum 


maximum 


maximum 


density 


squared 


Feedback gain for yaw damper 


sideslip angle 
yaw angle 
pitch angle 
roll angle 


distances metic 


distance in the 
distance in the 
flying tail 

aileron deflec- 


rudder deflec- 


Of sais 


divided by SF19 


Scale factor based on SF7 divided by SF1 and SF5 


Scale factor based on SF7 divided by SF1 and SF6 
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S33 
SF34 
Si35 
Sr 6 
SES 7 
SF38 
SI ee, 
SF40 


SFA 


SEACTOR 


SK 


SLOPE 


SN 


SSAn 


SX 


SXO 


SXN 


Sy 


SYN 


SZ 


Rate Jo reSh aco Sr 2 

Ratio Of SFI ters: 3 

SGake LAactOm basecevennhatlO7Of oh 5mEo Sr Ss 

Scale factor based on ratio of SF6 to SF9 

Scale factor based on ratio of SF6 to SF10 

Scale factor based on ratio of SF6 to SF1l 

scale factor based oOnjratio of SFO to ser 

Pheten@ @P =o 7. tor Sra 

Array Of scale £Lactors, made Up Of SES EhnriaseiZ 
Scale kactOnm Ck tie pia lacie. 


Variable indicating portion of each line to be 
removed from display in the window chop loop 


Slope of a display line, used in software window 
Negative of the scale factor 
Scalar equivalents of SA 


AGtual Cdistance Of alrecrake incl reocelen 
(amertia l ) 


initial Condl tien en ox 


Normalized distance of aircraft in X-direction 
(inertial) 


Actual distance of alrcraft im VaalrLecrelon 
(inertial) 


Normalized distance of aircraft in Y-direction 
ibiates ay eles a 


Detual cietance.oOf aircralt inw.z-d1recrion 
(inertial) 


inhi tral eCcOnalt lon on oc 
Final cutoff value on SZ (10,000 ft.) 


Normalized distance of aircraft in Z-direction 
(inertial) 


Mis Or Mie wale rat 
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TAU 


TLINE 


TM 


TRIG 


Bree GC 


vO 


VDOT 


VDOTO 


VINSCT 


VX 


VY 


VZ 


W 


x 


XEND 


XSTART 


YEND 


YSTART 


Time delay for the prediction and update 
Array of the earth grid reference points 
Thrust divided by mass of aircraft 
Array of sines of various angles 

Array of cosines of various angles 


A measure of roll angle in the X-direction used 
in generation of the "needle" 


A measure of roll angle in the Y-direction used 
in generation of the "needle" 


Velocity 

Initial condition on velocity V 

Time derivative of V 

VDOT of the iteration before 

Vertical intersection of a line used in window 
Welocity.in the X=-direction 

Velocity in sche: 2 -dlreetion 

Velocity in the Z4-direction 

Weight of the aircraft 

X value used in software window 

Array of X-coordinate ending points of display lines 


Array of X-coordinate starting points of display 
lines 


X-coordinate ending point for the horizon 

Temporary storage for X 

X-coordinate starting point for the horizon 

Y value used in software window 

Array of Y-coordinate ending points of display lines 


Array of Y-coordinate starting points of Gisplay 
lines 
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mE 


pre MP 


YTS 


ZEND 


ZOTART 


V Aad be oF 


ZTS 


Y-coordinate ending point for the horizon 


Temporary storage for Y-coordinate in software 
window 


Y-Coordinate Staremg poime for the horrzon 
Array of Z-coordinate ending points of display lines 


Array of Z-coordinate starting points of display 
lines 


Z-coordinate ending point for the horizon 


Z-coordinate starting point for the horizon 
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AEP EN DIT XD 


DATA DECK PREPARATION 


The following is a guide for use in preparing a data 
deck for use with the Ss ReMEOE: 

Cards 1 - 3 - These cards contain the inputed scale 
factors as defined in Appendix C. Each scale factor has a 


field width of ten and are ordered as follows: 


@erd) 1 — SFl, SF2, SF3, SF4, Sh5, She wGre) Gro 
@ard 2 - SF10, SF1l, SF12, SF13, SFl4, SFis5, selec 


Card 3 - SF18, SF19, SF30, the remaining columns are 
blank 


Sard 4 = This card contains the initial condi iens som 
the problem. Each number has a field width of ten and they 
are ordered as follows: initial X-coordinate (SXO) in feet, 
Mieeial altitude (SZO) in feet, initial velocity (VO) in 
ft/sec, initial angle of attack (Al10) in degrees, initial 
elevator angle (DITO) in degrees and the density of the air 


(eo) in Sita et: . The remaining portion of the card is 


plank. 
Cards 5 and 6 - These cards contain the aircraft con- 
stants. Each constant has a field width of ten and they are 


ordered as follows: 
Gard S = Weight (W) of airerait in 1b.7 wing span (eye 
ft, wing chord (C5 mine it. ,wine. aise dee 2) sern 
ft., moment of inertia about X-axis (RIXX) in 


slug-f£t*, moment Of 4nertiasabout the: Y-axis 
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(RI ai slug-ft“, moment of inertia about the 
Z-axis ({(REZZ) in slug-ft“, and the cross product 
of inertia (RIXZ) in slug-ft?. 


Cana 6 = Maximum thrust (T)-in lbs s. alti tudewmeutom: 


(SZF) in feet. The remaining columns are left 
blank. 
Cards 7 - 412 - These cards contain the aerodynamic 


coefficients dependent on both alpha and beta. The first 
card in each of a sequence of 135 three card series is as 
follows: in columns 1 and 2 the coefficient number, columns 
3 and 4 blank, column 5 the index on beta (IB), columns 6-10 
blank. The remaining portion of the card is filled with 
coefficients indexed on alpha (IA) in ten colum intervals. 
The second and third cards in the series contain the remain- 
ing coefficients in ten column intervals and indexed on 
alpha. Alpha is indexed 19 times for each beta index while 
beta is indexed 9 times for each of the fifteen coefficients. 
The range on alpha is 0 - 90 degrees with 5 degree intervals 


while the range on beta is -40 degrees to +40 degrees with 


iMmaegree intervals. The coefficients are as follows: 
Number Coefficient 
i Cy 
2 Co 
3 Co 
4 Cy 
rs) Cy. 
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Number Coe Ere term t 


Z 
7 E 
aie 
8 e 
; "Sit 
C 
10 “or ; 
C 
a are 
iat C 
Sr 
n2 C 
ae 
Nes C 
lea 
14 c 
Nea 
ES Cc 
*sit 
Cards 413 - 431 - These cards contain the aerodynamic 


coefficients dependent only on alpha. The first card in each 
of a sequence of 6 three card series is as follows: columns 
l and 2 the coefficient number, columns 3-10 blank, columns 
iieeo the coefficients indexed on alpha and in ten column 
intervals. The second and third cards contain the remaining 
coefficients again in ten column intervals and indexed on 


alpha. The coefficients are as follows: 


Number Coefficient 
16 (e 
Yp 
iy Cy 
p 
18 (& 
n 
Pp 
19 C 
YY 
20 Ci 
Yr 
2. C 
Nr 
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Cards 432-440 - These cards also contain aerodynamic 
Seetticients based on alpha only. The only difference in 
mie Setup of these cards in the field width is increased to 
eleven. As a eee tile the first card in the series is as 
follows: columns 1 and 2 the coefficient number, column 3 
mrank, COlumns 4-80 contain the coefficients indexed on 
eeepha and in eleven column intervals. The second and third 


cards contain the remaining coefficients indexed on alpha 


and in eleven column intervals. Phe coefiicients areas 
fhelLlows: 
‘Number Coefficient 
22 Cc 
Z 
q 
ZS ; Cc 
m 
S| 
24 Cy 
q 


Cards 440 -454 - These cards contain the actual coordi- 
nates for each line of the earth grid reference system. 
Bach card represents one line with its starting and ending 
coordinates. On each card the first eight columns are the 
starting X-coordinate, columns 9-16 the starting Y-coordi- 
nate, columns 17-24 the starting Z-coordinate, columns 25-32 
the homogeneous coordinate W, columns 33-40 the ending 
X-coordinate, columns 41-48 the ending Y-coordinate, columns 
49-56 the ending Z-coordinate, and columns 57-64 the homo- 
geneous coordinate W. The remainder of the card is left 


billank. 
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APPENDIX E 


OPERATING MANUAL 


Two tapes are available for use in execution of the 
digital program. One tape (labeled SPIN SIMULATOR ST) 
requires approximately 15 minutes to compile prior to execu- 
tion. The second tape (labeled SPIN SIMULATOR CD), which is 
a core dump of the first tape, requires only 10 seconds prior 
to execution. It .is recommended that the second tape be used 
due to its short compilation time unless a printed output of 
the program is desired. 

eae ot the required "ards, both control and “daranwoce 
located in the laboratory files under FIXED-BASE SPIN 
SIMULATOR. Also filed here is the master card version of 
the program and should not be used unless the tapes are 
destroyed. 

Detailed operating instructions for both the hybrid and 
graphics computers are available in the laboratory office. 
The following instructions should be followed in order to 
masure that the simulator is set up and executing properly. 

1. Load the tape selected on the tape drives. Select 

the proper mount (2 for SI, 3 for CD). 

2. Load and ready the card reader with control cards 

followed by the DATA. 

3. Ready the Line Printer. 

4. On the Digital Main Frame, Punch IDLF, RESET, RUN, 


CARDS. 
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Load the Analog Program board and the Analog Logic 
board into the computer (Boards number 7). 
On the Analog Keyboard punch KEYBOARD, LOCAL, POTSET. 
Im-the far right Gontrol box, only the following sen 
be lier]  1DCxXI or ADCS sland REALTIME Sa abl ao eines 
lights must be punched out. When complete punch 
DIGITAL (CMPTR. 
On the Analog the following options exist: 
DSO - UP = no printed output 
GUNTER = printed GutEpuE 
DS1 - UP - distance integrators disabled 
CENTER - distance integrators enabled 
DS2.-8UP = yaw damper ON 
CENTER I=" Vaw dampen Orr 
DOWN positions on these switches are mementary con- 
tact positions. 
At the Graphics terminal teletype to be used, load 
then Controlling progman by Eyoung 
PESEr ("GATDI , 104). ! 
Then execute the program by typing GATED! At this 
point nothing should be on the display screen. 
By this time the Digital computer should have 
compiled and loaded the main program and the follow- 
ing will be typed out on the digital consol teletype: 
INPUT AGT NUMBER 
The input light will then come on. To respond, type 
the AGT number (1 or 2) of the unit being used, 


followed by a carriage return. 


Se. 





1108 


Je 


ae 


3. 


ike 


Lae 


GC. 


ay 


BiG 


ye 


Zr. 


The Title should appear on the selected AGT. 
Tpeepatalwiill pe Tread in trem thereard. reader. 

If output is selected, the output will be printed 
on the line printer. 

The Analog Pots will be set by the Digital computer. 
The Address and Ratiometer readings should he 
echiangand . 

Set the cockpit in front of the AGP (wilieenr, 
reach to AGT—-2)2 The cockpit 1s “stored behancdmene 
Analog computer. 

The Instructions should flash on the screen. 

Punch the button on the Throttle plate and the 
display will appear. 

Punch the button on the control stick to fly. 1g 
during a run you wish to abort that run, pumen the 
button on the Throttle plate. 

At the completion of the run the Spin Results will 
be displayed, followed by a short delay in which to 
read them. 

At this time two options are offered. To fly again 
punch the: button cm the Control stack and thegeue. 
gram jumps back to the Instructions (#15). 

To receive the expanded program options, punch the 
Throttle button. Follow the displayed instructions 
and type your selection on the AGT teletype unit. 


This will result in one of three things happening. 
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i.. Logreal STOP to ehe proquam- 
2¢o JUMP tO, ene SnsStcuGtrons (5). 
3. Jump to the selection of the AGT number (#9). 
21. At the completion of the computer time, return all 
used equipment (cards, tape, boards, cockpit) to 


their storage locations. 
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eee NDT Xa 


SAMPLE OUTPUT 


The following is a sample output of the digital program 


which contains the following sections: 


1. Scale Factors and Aircraft Constants 

fee Original Alircratt Aerodynamic Cool tiecilenes 
3. Usable Aircraft Aerodynamic Coefficients 
Mee Pot settings for rtne Ane loq Comouter: 

5. Output of the Amplifiers and D/A Trunks 


6. Earth Grid Reference Lines 
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